Influence of the detector's temperature on the quantum Zeno eff'ect 
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In this paper we study the quantum Zeno effect using the irreversible model of the measurement. 
The detector is modeled as a harmonic oscillator interacting with the environment. The oscillator 
is subjected to the force, proportional to the energy of the measured system. We use the Lindblad- 
type master equation to model the interaction with the environment. The influence of the detector's 
temperature on the quantum Zeno effect is obtained. It is shown that the quantum Zeno effect 
becomes stronger (the jump probability decreases) when the detector's temperature increases. 
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I. INTRODUCTION 



The quantum Zeno effect is a consequence of the influence of the measurements on the evolution of a quantum 
system. In quantum mechanics the short-time behavior of the non-decay probability of unstable particle is not 
^ , exponential but quadratic jlj. This deviation from the exponential decay has been observed by Wilkinson et al. . In 
■ 1977, Misra and Sudarshan Q showed that this behavior when combined with the quantum theory of measurement, 
ly-j I based on the assumption of the collapse of the wave function, leaded to a very surprising conclusion: frequent 
1 observations slowed down the decay. An unstable particle would never decay when continuously observed. Misra 
and Sudarshan have called this effect the quantum Zeno paradox or effect. Later it was realized that the repeated 
measurements could not only slow the quantum dynamics, but the quanturnprocess may be accelerated by frequent 
measurements as well. This effect was called a quantum anti-Zeno effect |3;lH01- 
. The quantum Zeno effect has been experimentally proved in a repeatedly measured two- level system un derg oing 
--^ ■ Rabi oscillations. The interruption of Rabi oscillations has been at the focus of interest |1. 13. 1^ Eol ITU IT^ Il3l fT 3 . 
' Recently, the quantum Zeno effect has been considered for tunneling from a potential well into the continuum 
as well as for photoionization The quantum anti-Zeno effect has been obtained experimentally p^ . 

In the analysis of the quantum Zeno effect the hnite duration of the measurement becomes important and, therefore, 
the projection postulate is not sufficient to solve this problem. In Ref. ;6;] a simple model that allows to take into 
account the finite duration and finite accuracy of the measurement has been developed. However, this model does 
not take into account the irreversibility of the measurement process. 
IL^ ' The basic ideas of a quantum measurement process were theoretically expounded in Refs. IT^ l20ll2Tll2^ |23LI2^ 
on the assumption of environmentally induced decoherence or superselection. In this paper we extend the model, used 
in Ref. 6], including the interaction of the detector with the environment. Then it becomes possible to study the 
^ influence of other parameters of the detector on the evolution of the measured system, too. In this paper we analyze 
the influence of the detector's temperature on the measured system. 

To describe the decoherence and dissipation we use the the Lindblad-type master equation. Semigroup theories 
pioneered by Lindblad psl ] demonstrated that density-matrix positivity, translational invariance and approach to 
thermal equilibrium cannot be satisfled simultaneously. Under the assumption of Markovian dynamics and initial 
decoupling of system and bath, the semigroup approach adds dissipativc dynamics to the quantum master equations 
by means of the Lindblad dissipation operators. 

Recently the semigroup formalism has attracted much attention. Quantum computing is one of the fields in which 
quantum dissipation finds the most recent applications. In physical chemistry semigroup theories have been utilized to 
model dynamics of ultrafast predissociation in a condensed-phase or cluster environment '2(f\ , and electronic quenching 
due to the coupling of the adsorbate negative ion in resonance to the metal electrons in the desorption of neutral 
molecules on metal surfaces 27\. In nuclear physics, the semigroup formalism is applied to model giant resonances in 
the nuclear spectra above the neutron emission threshold p^ . 

We proceed as follows: In Sec. ^ we present the model of the measurement. The method of the solution is presented 
in Sec. Illll The measurement of the unperturbed system is considered in Sec. IIVI In Sec.[3we derived a formula for 
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the probability of the jump into another level during the measurement of the frequently measured perturbed system. 
Sec. I VII summarizes our findings. 



II. MODEL OF THE MEASUREMENT 



We consider a system that consists of two parts. The first part of the system has the discrete energy spectrum. 
The Hamiltonian of this part is Hq. The other part of the system is represented by Hamiltonian Hi. Hamiltonian 
Hi commutes with Hq. In a particular case the second part can be absent and Hi can be zero. The operator V{t) 
causes the jumps between different energy levels of Hq. Therefore, the full Hamiltonian of the system is of the form 
Hs — Hq + Hi + V{t). The example of such a system is an atom with the Hamiltonian Hq interacting with the 
electromagnetic field, represented by Hi, while the interaction between the atom and the field is V{t). 

We will measure in which eigenstate of the Hamiltonian Hq the system is. The measurement is performed by 
coupling the system with the detector. The full Hamiltonian of the system and the detector equals to 

H^Hs + Hd+Hi. (1) 

where Hjj is the Hamiltonian of the detector and Hj represents the interaction between the detector and the measured 
system, described by the Hamiltonian Hq. As the detector we use a harmonic oscillator with the Hamiltonian 

Hd = nn (^b^b + , (2) 

where b and w are the creation and anihillation operators, respectively. We choose the interaction operator Hj in the 
form 

Hi = XqHo, (3) 

where g = 6^ + 6 is the coordinate of the detector and the parameter A describes the strength of the interaction. This 
system — detector interaction is similar to that considered by von Neumann and in Refs. IsoIIsiIIb^ . Is^ Is^l . 
In order to obtain a sensible measurement, the parameter A must be large. 

The measurement begins at time moment to- At the beginning of the interaction with the detector, the detector's 
density matrix is poito). The detector initially is in the thermal equilibrium with the temperature T. Therefore, 



pnito) -PT= exp (^-^ j ^1 exp j j , (4) 

where fi = w b. The average excitation of the detector in thermal equilibrium with the temperature T is 

The full density matrix of the system and detector is /5(to) = ps{to) (g) poito) where ps(^o) is the density matrix of 
the system. 

The detector is interacting with the environment. The master equation for the density matrix of the system and 
the detector in the Lindblad form is (Ref. [25j ) 

^^l[H,p(t)]+Lnm], (6) 

where 

Loim] - E([^^'5(t), V;t] + [V„mV^]), (7) 

and are the Lindblad dissipation operators. We use the equation of a dissipative phase damped oscillator discussed 
in quantum optic [l^ . The Lindblad dissipation operators arc chosen as follows. 
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Then the equation (jBJ for the density matrix becomes 

+^(2atp(i)a - (n + l)p{t) - p{t){n + 1)) 

+^{2ap{t)a^ -hp{t)- p{t)n). (9) 
The approach to the thermal equihbrium is obtained when the parameters 7| and 7^ satisfy the condition [35| 

7T=7iexp(^-^j. (10) 

III. SOLUTION OF THE MASTER EQUATION 

For the solution of the equation Q we adopt the technique used in Ref. fs^. We introduce the quantum charac- 
teristic function [s^l 

x(C,r) = Tr{pe«^^-«*n- (11) 
The quantum characteristic function of the detector at the thermal equilibrium is 

XT{tn = eM~^Cn{T)). (12) 

We multiply the equation ijHJl by exp{—^*b) from the left and by exp(^6t) from the right and take the trace. When 
the interaction between the measured system and the detector is absent (i.e., A = 0), we obtain the equation 

d f d d 



7i 



t^e^x+r^xj. (13) 

We will search the solution of Eq. in the form 

Xitn = exp (^C.Mm'i-C)"^ , (14) 

where Cj^k are the coefficients to be determined. Substituting Eq. (|14|l into Eq. H13|l we obtain the set of equations 
for the coefficients Cj^k 

+ ^(7T - 7i)(j + fc)Cj,fe(i) + 7t'5jm4,i- (15) 

The solution of Eq. is 

Ci,i(t) Ci,i(0)e-(''i-''T)t + ^(T) (1 - e-^'^i-'^T)*^ ^ (16) 

Qfe(t) = Cj-fe(0)e'"(^-'=)*e-*(^-*=)'*-5(Ti-'^T)0+fc)*, j^l,/c^l (17) 

From the solution we see that the function x approaches the function at the equilibrium xt as the time t grows. The 
detector's density matrix p}), correspondingly, tends to the pr- 
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IV. MEASUREMENT OF THE UNPERTURBED SYSTEM 

At first, we will consider the case when the perturbation is absent, i.e., V(t) = 0. Since the Hamiltonian of the 

measured system does not depend on t we will omit the parameter in this section. We can choose the basis \na) 
common for the operators Hq and Hi, 

Ho\na) = E^\na), (18) 

Hi\na) = Ei{n,a)\na), (19) 

where n numbers the eigenvalues of the Hamiltonian Hq and a represents the remaining quantum numbers. 
We introduce the density matrix pm.n = ^^'^ characteristic function 

XmAi,C\t) = Tr{p„,„(t)e«^'e-«*^}. (20) 
From Eq. ((SJ we obtain the equation for the density matrix pm^n 

Q 

-Q^Pm.n = —i^mnPm.n — i^{oJmqPm,n — Pm,nq^n) + L]j[pm.,n], (21) 

where 

= (22) 

Wmn — ~ ^n- (23) 

Equation (|21|l may be solved similarly as in Sec. lIIII When the detector is initially at equilibrium then Xmn{£,, 0) = 
Xmn(O) cxp(— ^^*n(T)). As in Sec. IIIII we take the characteristic function of the form (fTHl and obtain the equations 
for the coefficients Cj^k 



dt 
dt 
dt 



-lUmn (1 + A(Ci,o + Co,l)) , (24) 

{in - 7cff)Ci,o + «A(tJ„ - ujmnCis), (25) 

-{in + 7cff)Co,i - iX{oj,n + w,„„Ci,i), (26) 

(7T -7i)C^i,i +7T (27) 



dt 

with the initial conditions Co,o(0) = 0, Ci,o(0) = 0, Co,i(0) = 0, Ci,i(0) = n{T). Here 



7eff = ^(7 + 7i -7t)- (28) 



The solutions of Eqs. (|2IIl-l|2ZI) are 

Cia(t) = n{T), (29) 

C, ,{t) = ^^ ^n-u;^nn{T) / ^ _ ^(,a-y.«)A (30) 

7cff ~ iil \ / 
Coi(t) = _^^ ^-^ +^^^nn{T) (31) 
7off + \ J 

CoAt) - + A^^„,„ "" ~ "-.^^^^ (t + - l)) 

7off - III V 7cfr - V A 

-A^c.^„ ""+'^-^^^^ / 1^ _ y (32) 

7off + 252 V 7cff + ^ ^ / 

Using Eqs. (|14|l and (|32|l we find that the non-diagonal elements of the density matrix of the measured system become 
small as the time t grows. This represents the decoherence induced by the measurement. The diagonal elements of 
the density matrix do not change. 
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V. MEASUREMENT OF THE PERTURBED SYSTEM 

The operator V{t) represents the perturbation of the unperturbed Hamihonian Hq + Hi. We wiU take into account 
the influence of the operator V by the perturbation method, assuming that the strength of the interaction A between 
the system and detector is large. 

The density matrix at time t is related to the initial density matrix by the equation p{t) — S(t)p[0). The superop- 
erator S obeys the equation 

§iS = LS, (33) 
where the Liouvillian superoperator L is defined by the equation 

Lp^^[H,p]+Ld[p]. (34) 
in 

Here H and Lo are defined by Eqs. and {Tj), respectively. We can write L = Lq + Ly, where Ly is a small 
perturbation, defined by the equation 

LvP=Uv,p]. (35) 
tn 

We expand the superoperator S into powers of V 

S = ^(0) + 5(1) + 5(2) + . . . . 

Then from Eq. (ESJ it follows 

|5(°) = LoS("\ (36) 
ot 

|-5(*) = Lo5(*) +Lv5(*-i). (37) 
ot 

The formal solutions of Eqs. H36() and H37|l are 

5(°) = e^°* (38) 

and 

5(')= /*dii5(°)(t-ti)Ly5(''-i)(ti). (39) 
Jq 

In the second-order approximation we have 

S{t) = 5(")(t)+ / dii5(°)(i-ti)iy5(°)(ti) 



Jo 

+ f dh r dt2S^''\t-h)LvS^°\ti-t2)LvS^'>\h). (40) 
Jo Jo 

Let the initial density matrix of the system and detector is 

p{Q) ^\ia){ia\(^ pD, (41) 

where p'o is the density matrix of the detector. The probability of the jump from the level \ia) to the level |/ai) 
during the measurement is 

W{ia ^ /ai,t) = Tr{|/ai)(/ai|p(i)}. (42) 
The unperturbed evolution does not change the energy of the measured system, therefore, we can write as 

5(*')(t) [\ma){na'\®pD] = \ma){na'\® s':^l^^^,{t)pD- (43) 
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Equation H43|l defines a new superoperator S'^^ acting only on the density matrix of the detector. The indices 

ma and na' in 5*^^^ denote the states of the measured system. From Eq. (|^ it follows that the superoperator 

Sma,ma with cqual indices does not change the trace of the density matrix p/j, since the trace of the full density 
matrix of the measured system and the detector must remain unchanged during the evolution. 

We assume that diagonal matrix elements of the perturbation operator V are zeros. Inserting the expression 
p{t) = S{t)p{0) into Eq. 142|l and using equation l|40|l for the superoperator S{t), we obtain the jump probability 



W{ia^fai,t) = 1^ I dh I dt2TT{\fai){fai\ 



-S^"\t - h)[S^"\h - t2)[S^"Ht2)m]V]v) } . (44) 

From Eqs. H43|) and (|41|) it follows that the last two terms in Eq. (|44|l contain the scalar product {fai\ia). Since the 
states l/ai) and \ia) arc orthogonal, the last two terms in Eq. H44(l are zeros. Therefore, the jump probability is 

W{ta-^fai,t) ^ 1^ f dh / 'rft2Tr{|/ai)(/ai| 
n Jo Jq 

+ s^°Ht-ti)[s^°Hti~t2)vs^°Ht2)pmv)}. (45) 

From Eq. H45(l . using expression for the initial density matrix of the system and the detector (|41(l and equation H43|l . 
we have 

n Jo Jo ^ 

X - t2) + Sfl^^jh - t2)) S^Zjt^)PD} ■ (46) 

The superoperator Spl_^ {t — ti) preserves the trace of the detector's density matrix, therefore, the jump probability 
equals to 

iy(za^/ai,i) = ^iV^cjaA^ J^h j\t2Tv[(^s'll]f^^{h~t2) 

+ sZ,c.iti-t2))slZMPD}- (47) 
Defining a new characteristic function similarly as in Eq. (|2U|I 

X..Jo.(^,r;ii,i2) = Tr{e«^^-«*^^15/„,(ti-<2)5l°),Ji2)p^ (48) 
the jump probability (|47|l can be expressed as 

W{ia^ fai,t) ^ ^\V,aja,\^ [ dti [ (x»a,/ai (0, 0; ti, tz) + x/ai,*a(0, 0; ti, t2)) • (49) 







The detector initially (at t = 0) is in the thermal equilibrium with the temperature T, pu = pr (Eq- O)- The 
initial characteristic function is Xia./ai (■fj C*; 0, 0) = exp (— ^^*n(T)). Using the results of the Sec. IIVI fEas. ((TH) and 
()29 |l -l|32 |) ), we obtain the characteristic function of the density matrix at time ^2 

X^c.J»AC,C■,t2,t2) - Tv{e^''''e-i''sl°],^{t2)pD} 

VToff - «i2 V J J 

X exp ( J^I^ ( 1 _ e-'(^^+^"«)*0 - ^Cfi{T)) . (50) 
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Taking the function Xiajai (Ci C*! ^2, ^2) from Eq. (|50|l as the initial characteristic function and preceding further as 
in Sec. IIVI we have the value of the characteristic function, defined by Eq. (|48|) with the parameters ^ = ^* = 



X^aJa^iO,0■,tl,t2) = CXp iuJiaJaiih - 

7off - «ii V 7cff - «S2 V / 

7cff + iQ \ 7cff + lii V 

I -^^^''Z^' ^ _ g(iO-7„tt)t2'j ^ _ g(*n-7off)(«l-t2) 

(7off - V J \ 



A UJj fUJ. 
ilcS 



(51) 



Here 

1 



WiQjQi = t^i/ + - - £'i(/, ai)) . (52) 

A. Approximations 

When the dissipation is fast, i.e., the dissipation time is much less than the period of the oscillator, we have f2 <§; 7eff . 
Then 

Xia,/ai(0, 0;ii,i2) = exp (-iwiQjQi(ii - ^2)) 

X exp ( -(1 + 2n{T))^^ {h -t2 + — f e-^=«(*^-*^) " M ■ 
\^ 7cff V 7cff ^ ^ J j 

The probability of the jump from the level \ia) to the level \fai) during the measurement according to Eq. H49|l is 

/ (l + 2n(r))AV^ f 1 \\ 

X exp ^ u H (e - 1) • (53) 

\ 7efr V 7cff J ) 

We introduce the function 

Ht) fau'^c = \V^c,fc., I' exp (^^[Eiif, ai) - Ei{t, a)]?j (54) 
and the Fourier transformation of <i>(t)/Qi,j;a 

1 f°° 

G{uj)fai.ia = TT dt<P{t)fai.iaexp{-lUjt). (55) 

Then we can rewrite Eq. (|53|l in the form 

W{ia ^ fai,t) ^ —Y dijjG{uj) fai^iaP{^)if , (56) 

where 



P{uj)if — —Re [ du(l ^ exp {i{Ljj — LUif)u) 



/ (1 + 2n{T))\^Ljff /I 
X exp — ^ u + _(e-T"«" - 1) . (57) 

V 7eff V 7eff 
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The equation 1)56(1 is of the form, obtained by Kofman and Kurizki |5'|, assuming the ideal instantaneous projections. 
The function P{uj)if is the measurement-modified shape of the spectral line (Refs. (5i. i6i. i38i| ) . Here we have shown 
that Eq. ((56(1 can be derived from more realistic model as well. The assumption that dissipation is fast, Vl <^ 7cff is 
crucial. Without this assumption the jump probability cannot have the form of Eq. 1(56(1 . since then Xia.fai (0, 0; ti, t2) 
depends not only on the difference ti — t2 but also on t2- 

When A is big then to the integral in Eq. 1(53(1 contribute only small values of u and we can expand the exponent 
exp(— 7cffM) into Taylor series keeping the first three terms only. We obtain the jump rate 

R{ia fai) « -^iViajaA^ Re J du exp (^ujfai,iaU - ^(1 + 2n(T))A^w5^u^ j 

or 

""^^^ ^"^^ " -WkJV 2(i + 2f^(r)) 

The obtained decay rate inversely proportional to the measurement strength A. The measurement strength appears 
in the equations multiplied by i/l + 2n{T), therefore, the effect of the measurement increases as the temperature of 
the detector grows. 



VI. CONCLUSIONS 



We analyze the quantum Zeno effect using the irreversible model of the measurement. The detector is modeled 
as a harmonic oscillator, initially being at the thermal equilibrium. The interaction of the detector with the system 
is modeled similarly as in Ref. (6||. The Lindblad-type master equation for the detectors density matrix is solved 
analytically. An equation for the probability of the jump between measured system's states during the measurement, 
similar to that of Refs. is obtained ((56(1 . From the used model it follows that the increase of the detector's 

temperature leads to the enhancement of the quantum Zeno or quantum anti-Zeno effects. 
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